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PREFACE 


This  report  describes  work  done  in  the  summer  of  1983  by  Dr.  B.  D. 
Slvazlian,  Department  of  Industrial  and  Systems  Engineering,  the  University  of 
Florida,  Gainesville,  Florida  32611  under  Contract  No.  F08635-83-C-0202,  with 
the  Air  Force  Armament  Laboratory  (AFATL),  Armament  Division,  Eglin  Air  Force 
Base,  Florida  32542.  The  program  manager  was  Mr  Daniel  A.  Mclnnis,  (DLYW) . 

The  work  was  initiated  under  a  1982  USAF-SCEEE  Summer  Faculty  Research 
program  sponsored  by  the  Air  Force  Office  of  Scientific  Research  conducted  by 
the  Southeastern  Center  for  Electrical  Engineering  Education  under  Contract 
No.  F49620-82-C-0035. 

This  work  addresses  itself  to  the  problem  of  computing  the  uncertainty 
associated  with  the  probability  of  kill  P^  due  to  fragmentation  in  the 
presence  of  aiming  error  and  in  the  absence  of  blast.  Let  P^  be  the 
probability  of  kill  due  to  fragmentation  in  the  absence  of  aiming  error. 

Assume  that  points  on  the  ground  surface  are  referenced  relative  to  a  system 
of  Cartesian  coordinates  where  the  x-axis  is  pointed  in  the  direction  of  range 
and  the  y-axis  is  pointed  in  the  direction  of  deflection.  The  origin  0  is 
arbitrarily  selected.  For  a  point  target  located  at  (x,y)  assuming  that  the 
weapon  bursts  at  a  point  (u,v),  the  probability  of  kill  due  to  fragmentation 


The  aiming  errors  in  the  direction  of  range  and  deflection  are  assumed  to  be 

unbiased  centered  at  (x,y)  independent  of  each  other  having  a  Guassian 

distribution  with  respective  standard  deviation  a  and  a  . 

x  y 

It  is  shown  that  can  be  expressed  as  a  mathematical  function  of  the 

four  parameters  Rx,  R  ,  o  and  a  .  Moreover,  under  the  assumption  that  Rx, 

y  x  y 

R  ,  a  and  a  are  not  known  with  certainty  but  are  estimates,  explicit 
y  x  y 

expressions  are  obtained  for  E[P|<f]  and  VarfP^]. 

The  author  has  benefited  from  helpful  discussions  with  several  people. 
Particular  thanks  are  due  to  Mr  Jerry  P.  Bass,  Mr  Daniel  A.  Mclnnis  and 
Mr  Charles  A.  Reynolds,  all  from  DLYW,  who  have  read  the  report  and  who  have 
contributed  to  it  through  helpful  comments. 

The  report  is  the  third  of  a  series  of  four  reports  dealing  with  the 
uncertainty  associated  with  various  weapon  effectiveness  indices,  and  details 
methodologies  and  techniques  used  in  computing  such  uncertainties  in  the 
presence  of  error  in  the  input  parameters. 

The  Public  Affairs  Office  has  reviewed  this  report,  and  it  is  releasable 
to  the  National  Technical  Information  Service  (NTIS),  where  it  will  be 
available  to  the  general  public,  including  foreign  nationals. 

This  technical  report  has  been  reviewed  and  is  approved  for  publication. 
FOR  THE  COMMANDER 

MILTON  D.  KINGCAID,  (gblonel,  USAF 

Chief,  Analysis  and  Strategic  Defense  Division 
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SECTION  I 


INTRODUCTION 

This  report  considers  the  problem  of  estimating  the  probability  of  kill 
due  to  fragmentation,  Pkf,  in  the  presence  of  aiming  error.  Prior  to  solving 
the  estimation  problem,  it  is  necessary  to  derive  a  mathematical  expression 
for  Pkf  in  order  to  apply  the  usual  statistical  techniques  to  arrive  at 
confidence  intervals  for  Pkf.  In  Section  II,  the  derivation  of  Pkf  is 
endeavored  based  on  several  explicitly  stated  assumptions.  In  Section  III, 
the  estimation  of  E[Pkf]  and  Var[Pk^]  is  carried  out  when  using  the  subjective 
estimation  procedure.  In  Section  IV  the  estimation  technique  using  Taylor's 
series  is  applied  to  obtain  E[Pkf]  and  Var[Pkf].  Finally,  Section  V  provides 
some  conclusive  remarks. 

A  comparison  of  Section  III  with  Section  IV  will  show  that  the  subjective 
estimation  procedure  involves  complex  mathematical  expressions  which  result  in 
cumbersome  calculations.  Furthermore,  the  procedure  assumes  that  the 
estimates  are  independent  random  variables.  This  assumption  is  not  necessary 
when  using  the  Taylor's  series  estimation  technique.  In  addition,  the  latter 
allows  one  to  segregate  the  contribution  of  each  variance  component  to  the 
total  Var[Pkf]. 
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SECTION  II 


MATHEMATICAL  MOOEL  FOR  Pfcf 

1 .  Background 

In  this  section,  a  mathematical  model  for  the  probability  of  kill  due  to 
fragmentation,  p^,  in  the  presence  of  aiming  error  is  developed.  The  basic 
situation  that  one  is  facing  consists  of  the  following. 

A  weapon  whose  main  effect  is  kill  due  to  fragmentation  is  delivered  from 
.it  to  a  target  point  located  on  the  ground  surface.  The  target's  position  is 
stationary.  The  weapon  may  not  directly  hit  the  target  due  to  the  presence  of 
aiming  errors.  These  errors  are  assumed  to  be  unbiased;  that  is,  centered  at 
the  location  of  the  target.  The  target  may  or  may  not  be  killed  by  the  effect 
of  fragmentation ;  thus,  the  target  is  killed  with  a  given  probability  level. 
The  probability  of  kill  due  to  fragmentation  is  related  to  the  distance 
between  weapon  and  target  by  a  well-defined  mathematical  function.  In 
addition,  the  aiming  error  is  not  known  precisely  but  is  expressed  by  a 
probability  density  function  which  provides  a  mathematical  formula  for 
computing  the  probability  that  the  weapon  will  impact  in  an  interval  du  dv 
:  ’  **  s**  fn  a  point  (u.vl  on  the  ground  surface. 

The  technique  that  will  be  used  to  compute  the  expression  for  the 
probability  of  kill  due  to  fragmentation,  Pkf,  in  the  presence  of  aiming  error 
is  based  on  the  laws  of  conditional  probability.  Ultimatel  y.  Pkf  is  not  going 
to  depend  on  the  position  of  the  target  if  the  aiming  error  is  unbiased,  and 
if  weapon  delivery  can  theoret ical ly  result  in  a  point  of  impact  which  can  be 
3'Vwhore  on  the  ground  surface.  On  the  other  hand,  will  depend  on: 

a.  the  parameters  specifying  the  functional  form  relating  probability  of 
kill  to  di stance ; 

b.  the  statistical  parameters  of  aiming  error  distribution. 

? 


2.  Assumptions 

The  following  assumptions  pertaining  to  this  situation  will  be  made: 

a.  Each  of  the  target  and  weapon  is  idealized  as  a  point,  and  the  weapon 
is  aimed  at  the  target. 

b.  The  direction  of  the  weapon  delivery  range  and  deflection  are, 
respectively,  parallel  to  the  (x-y)  coordinate  system  on  the  ground  plane. 
Since  the  coordinate  system  can  be  arbitrarily  selected,  there  is  no  loss  in 
generality  •  '  making  this  specific  assumption.  The  position  of  the  target  has 
coordinates  (x,y). 

c.  The  aiming  error  (distance)  in  each  of  the  x  and  y  directions  are 

independently  and  normally  distributed  with  respective  means  x  and  y  and 

standard  deviation  a  and  a.  Let  (du,  dv)  be  the  infinitesimal  rectangle 

*  y 

close  to  the  point  (u,v)  at  which  the  weapon  impacts,  and  define  the  random 
variables  U  and  V  which  measure,  respectively,  the  distances  between  the 
target  point  and  the  weapon  impact  point  along  the  abcissa  and  the  ordinate. 
Then,  the  probability  that  the  weapon  will  impact  in  the  rectangle  du  dv  is 

fU,V  (u-x’  du  dv  =  exP[-  i72L]  * 

x  y  ?.o 

J  x 

2 

exp[-  ^-v"^  ■]  du  dv  .  (1) 

2oy 

d.  The  probability  of  kill  due  to  fragmentation  at  a  point  (x,y)  given 
that  the  weapon  impacts  at  (u,v)  is  given  by  the  Carleton  damage  function 

2  2 

%<*-“.  y-v)  5  ex?l-  (tt)  ■  (*ir)  1- 

x  y 
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f*iis  is  sometimes  known  as  the  elliptical  damage  function  involving  the  two 
parameters  R^  and  R  These  parameters  are  often  identified  as  the  weapon 
radii:  Rx  is  known  as  the  range  weapon  radius,  and  R^  is  known  as  the 
deflection  weapon  radius. 

e.  Blast  effect  is  neglected.  This  implies  that  the  target  is  not  blast 
sensitive.  Or,  if  weapon  blast  exists,  its  effect  on  the  target  is  negligi¬ 
ble,  hence  not  resulting  in  a  kill. 

f  .  In  general  ,  the  point  in  space  from  which  weapon  is  delivered  is 
nor,  st  at  i  on  ary ,  the  weapon  is  subject  to  ballistic  errors,  etc.  It  shall  be 
assumed  that  all  these  factors  combine  into  a  single  source  of  error  which  is 
incorporated  in  the  aiming  error. 

g.  The  fragmentation  does  not  contribute  to  the  aiming  error. 


1 .  The  Mathematical  Model 

The  probability  of  kill  at  (x,y)  due  to  fragmentation  is: 

00  oo 

P|<f  =  /  /  [Probability  of  kill  at  (x,y)  |  weapon  impacts  at  (u,v)] 

—  OO  -OO 

[Probability  that  the  weapon  impacts  between  (u,v)  and 
(u+du,  v+du)1 


o 


kf 


x>  oo 

/  / 

—  OO  — OO 


o  cr» 

/ 

—  c*>  —  oo 


1 

.3  a 

x  y 


Pk(x-u,  y-v)  •  f(JjV(u-x,  v-y) 


expt 


La --*!-] 


du  dv 


du  dv  . 


Maxing  the  change  in  variables  w  =  u-x  and  z  =  v-y  yields 


hi  ■  /  /  “Pi-  •  ST-"11!-  [7T*^7)I  dW  dZ  ' 

-»  -»  Rx  Ry  x  y  2  x  2oy 


This  integral  may  be  expressed  as  the  product  of  two  single  integrals  as  follows: 


pkf  =  J 

/  2n  0  -“> 

X 


1 - /  exp[-  (— \  +  —-?)  W^J  dw 


f-?  — 1> 
R  2< r 

X  X 


.  - 1 - /  exp[  -  (-1  +  -Lj  Z2j  dz  . 

-»  Ry  2  oy 


The  first  integral  is  computed: 


L  =  - - - /  exp[-  (-1  +  -U)  w2i  dw 

/-2Fax  -»  R2  2o2 


Ti'  /7  to* 

/  x  x 


I  e 

x  /~?T  a  -- 
x 


~T  +  “ ~7 
<  2°x 


Since  J  e  2  de  =  /2jt 


It  fol lows  that 


o  n  rr~7~r 

x  ~Z  — Z 

J  S  2°x 


Rx  +  2°x 


Similarly,  one  obtains  for  the  s econrl  integral  in  (3) 


ly  = 


/TtT 


/  ’  exp[  -  (-£  +  -L)  z 
R  2a 


y 


Substituting  (S'  and  (6)  in  (3)  yields 


pkf 


I  I 


R 

_ x 

- X - X- 

/  R^  +  ?o^ 


R 

y _ 

X-  * 

✓  R  +  2ac 


dz 


SECTION  III 


ESTIMATION  OF  E[Pkf]  AND  Var[Pkf]  USING 
THE  SUBJECTIVE  ESTIMATION  PROCEDURE 


1.  Background 

Relation  (7)  shows  that  the  mathematical  expression  for  Pk^  is  a  function 

of  the  four  parameters  R  ,  R  a  ,  and  a  .  In  the  subjective  estimation 

*  y  x  y 

procedure,  it  is  assumed  that  the  uncertainty  level  of  each  of  the  four  input 

parameters  is  provided  as  a  subjective  information.  A  lower  and  upper  bound 

value  for  each  parameter  is  obtained.  These  may  be  determined,  for  example, 

through  a  subjective  procedure  in  which  individuals  are  requested  to  provide  a 

lower  and  upper  bound  values  on  R  ,  R  ,  a  ,  and  a  based  on  their  judgement 

y  x  y 

and  their  experience.  The  value  of  a  particular  parameter  is  assumed  to  take 
on  equally  likely  values  between  its  two  extreme  points.  This  is  equivalent 
to  assuming  that  each  parameter  is  a  random  variable  uniformly  distributed 
over  its  range  of  values.  Further,  the  parameters  are  assumed  to  be  mutually 
independent  random  variables.  With  this  statistical  information,  the  evalua¬ 
tion  of  E[Pkf]  and  Var[Pkf]  are  reduced  to  the  computation  of  a  set  of 
definite  integrals. 

We  now  assume  that  the  four  parameters  Rx,  Ry,  o  ,  and  a  are  estimated 

j  x  y 

independently  and  uniformly  distributed  over  the  following  ranges: 


o 

X 


1 

1 


<  R  <  R 

X  X, 


<  o  <  a 
x  x. 


2.  Estimation  of  E[ P k 


<  R 

y 

<  a 


y 


<  R 

y 


<  a 

y 


2 

2 


From  relation  (7),  it  is  clear  from  the  stated  assumptions  that  I  and  I 

*  y 

are  independently  distributed.  Hence, 
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E[Pkfl  =  E[lxl  PL  i  y  J  •  (8 

Thus,  the  problem  reduces  to  finding  E[I  ]  and  E[I  ].  To  simplify  notations, 

x  y 


/  R 2+2a2 

wh ere  R  and  a  are  independently  and  uniformly  distributed  over  the  respective 
ranges  R^  <  R  <  R?  and  a  <  a  <  a^.  If  E[  II  is  calculated,  then  EC  I  X"1  and 
Efl  ]  can  be  immediately  determined.  Now 


fii  /iv 


dR  a  /  d(/  R  ^  2a‘ 


R^  +  2a2  -  /  R2  +  2a2  . 


Substituting  (11)  in  (10)  results  in 


-Tp— s-fr - Y  I  (/ Ro  +  2  2  -  J  r\  +  2a2)  do 

(R2-Rl)(°2'°l)  o  2  1 


R1  7 


i  [»  *  (4i 

/  o 


Using  (16)  in  (15)  one  obtains 


1 


ECU  = - 

y  ^T(RU  -RU)K  )  J2 


y2  yf  y2 


L0y,  (Y22“Y12^  ■  °y1  ^Y21'Y11 


y2  ,  y2  ^1+Y22^ 
+  -ln^Ti^y 


R 


yl  ,  y2  ^ 1 +Y 1 2 ^ 


1  *—  i  <- 

I  n  — —  tttu — rj  • 


T  '"^TT^7 


3.  Estimation  for  Var[P^] 

From  relation  (7)  one  obtains 

Var[Pkf]  =  E[I2]  •  E[I2]  -  |E[IX]  •  E[Iy]}2  . 

o 

Since  E[IX]  and  E[Iy]  have  been  obtained,  it  suffices  to  find  EC  1^3  and 
L [ 1 2 ] .  Expression  (9)  is  referred  to  and  one  has 


Thus: 


1*2  <*2  2 

ED?]  -  in  u  Ho  ,|f  ! 

r<2  Vl°2  ai>  R1  0l  R2+2o2 


(19) 


Now 


r°2  R2  .  R2  f°2  do 

J  -9 - 7  da  =  -r-  I 


~2 - 7  uo  ~  "T-  ■»  “5  7 

°i  R  +2°  °1  °  +  (JL) 

n 


R  1  a 

7“  33  1  HT 

✓~2  /~2 


'1 


R  -1  °2  R 

- 1 3n  —  ~ ~p  "  •"  —  — 

n  K  rr 


tan 


■  1  'X™ 


JL  cot'1  -JL 

n 


—  cot"^ 


n 


(20) 


Substituting  (20)  in  (19)  yields 


cr  T  2n  1  f  [  R  *.“1  R  R  “  1  R  ]  jq  /o,\ 

E[I  ]  =  — n  -yra — ~ — r  /  L—~  c°t  — ~  -  — —  COt  J  dR  .  (21) 

(K2  K1M°2  °1}  Rx  ozT7  n.  o 


Making  use  of  the  following  result 


fx  (cot-1  ax)  dx  =  cot'1  ax  +  i-  +  const, 

2a2  23 


2  2 

1+aV  _„-l 


expression  (21)  becomes 


13 


■ .  .J,  .  II  ■  vu-  vj  i  '.-l  \t  '7»  '.V  1 i.  I  .V-.  rTTWT-WT-T-  v- T- 


E[rl  = 


i 


<VV(vV  /-r  2-L 

2o2 


1  +  R2 

[ - ^ -  cot’1  R 


Or,  * '2 


1^«2 

201  -1 
-  cot  1 


1 


a2  n 


2  a 


/T  r2+2o2 

TR'rTi)(°2-°i)  tdn 


R2+2o?  .  a.  AT  a.R 

- r-1  tan"1  - L -  1 


R 


2/2 


This  expression  finally  reduces  to 


E[I2] 


2  2 
Ro+2a, 


at  r  < 

(R2-Rl)(a2-a1V  I  4 


tan 


R 

R 

oj  AT 

2  — - - 

0^  AT 

0?  /  2 

R 

T 

2  /"T 

R2 

]  • 

R1 

a2  AT 

rv 

2  2 
R2+2°  1  f 

o,  AT 


tan 


R2+2o2  ,  o9  AT  R?+2a?  .  o.  AT 

tan"1  A - +  J-rJ-  tan"1  A - ]  ♦  | 


(22) 


V1 


V1 


2  2 

The  expressions  for  E[I  ]  and  E[I  ]  are  immediately  obtained  as 

x  y 


AT 


2  1 

E[Ix]  =  2  +  ( R v  -Rv  )(av  -a  ) 


R2  +2o2  0  AT 

_ T  t  -1  X2 

[  A  t3n  D 


x2  x1  X?  Xj 


R2  +2o 2  0  /  2 

X2  X1  t  -1  X1 

A  1 3  n  n 


R2  +2o2  o  /  2 

X1 _ X_2  t  -1  X2 

a  l  3  n  n 


R2  +2o2  o  /  2 

♦  tan'1  ■ 


R 


(23) 


kl 
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•— 


?  2 
K  +2cr  . 


o  /  2 


Ff  1^1  -  —  +  ^  2 

Et  V  "2  (R  -R  )  (o  -a  ' 

y2  yl  y2  yl 


-1  y2 


R2  +2o2 

y2  yj 


.-l  yi 


R2  +2o2 


VI 

\ 


R2  +2o2 


a  /  2 


.-1  yl 


4.  Example 

A  weapon  whose  main  effect  is  kill  due  to  fragmentation  is  aimed  at  a 
particular  target  from  a  flying  aircraft.  Following  its  release,  it  is  estimated 
that  the  weapon  impact  angle  will  he  60°.  The  standard  deviation  of  the  aiming 
error  in  range  is  150  +_  50  ft  (to  be  interpreted  as  equally  likely  between  100 
and  200  ft).  The  standard  deviation  of  the  aiming  error  in  deflection  is  80  +_  10 
ft.  For  the  given  impact  angle  and  weapon/target  situation,  it  is  estimated  that 
the  range  weapon  radius  and  the  deflection  weapon  radius  are,  respectively ,  85  _+ 

5  ft  and  170+^8  ft.  It  is  required  to  determine  the  probability  of  kill  Pkf  of 
the  weapon  and  to  provide  a  two-standard  deviation  confidence  interval  for  Pk^. 
Using  the  formally  introduced  notations  one  has 


R  =  80  ft,  R  =  90  ft,  R  =  162  ft,  R,  =  178  ft, 
1  2  yl  y2 


o  =  100  ft,  a  =  200  ft,  o  =  70  ft,  and  a  =  90  ft, 
1  X2  yl  y2 


From  equation  (13)  the  values  of  Xi i  are  first  derived 


'll 


i  .  A  ♦  J22L 


2o 


2(100)  ‘ 


12 


f  \  / - 2 

i  *  -  /7*Zmi . 


2a 


2(200)' 


X21  =  7  1  + 


^  =  /i  +Ji2L 


2a 


2(100) 


22 


i  *  J^./i 


7" 


2a2  2(200) 2 

X2 


From  equation  (16)  the  values  of  are  next 


R 


T 


-/ 1 .  .  A  * ml 

2o2  2(7  0) 2 

yl 


/.  y 

Y12  =  ^  1  +  — ? 

2a  ^ 


‘  .  47  ml . 


2(90) 


R 


Y.. 
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. .  2* .  47ml 

2o^  2(7  0; 

yl 


Y22 


1  + 


/um!-  i 


2o‘ 


2(90) 


1.148,912,6 


1.039,230,5 


1  .185,326,9 


.049  ,404,6. 


derived 


.917,800,6 


.618,641,4 


.057  ,440,5 


.719  ,244  ,7. 
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— c  -n - «  -T --V  -  -V  -  V  -  v*  V  »-.  *  l  ■  V  »  ■-  'T'T'  -  ’ 


L'.V 
. ' 

L'  *  ’ 


ft 


Referring  to  (14),  E[IX]  is  computed  using  the  numerical  values  of 
Xj j  (i.j  =  1,2), 


ECU  - 


1 


r?  (90-80)  (200-100) 


l  (200)  ( 1 .049  ,404,6  -  1.039,230,5) 


-  (lOOp  (1.185,326,9  -  1.148,912,5) 


x  (90p  .  200  1  +  1.049,404,6 

~T~  100  1  +  1.185,326,9 


,  (30)  200  1+1.039,230,5, 

+  — jr-  in  THD  I+T;T4’8,gr2,'5:J 

1  [(40, 000)  (.010, 174,1) -(10, 000)  (.036, 414,4) 


E[U  = - 

X  /“Z  (10)  (100) 


+  (2I°£j  (.628,931,1)  -  (ii£22)  (.640,757,8)] 

=  .381,530,8.  (25) 

Referring  to  (17),  E[Iy]  is  computed  using  the  numerical  values  of  Y^j 

(i,j  -  1,2). 


ECU  = 


y  r~Z  (178-162)  (90-70) 


[ (90;  ( 1 .719  ,244,7  -  1.618,641,4) 


(7  0) ^  (2.057,440,5  -  1.978,006,0) 


♦  iizailn 


90  1  +  1.719,244,7 

7U  1  +  2.057,440,5 


(162)  i  _  90  1  +  1.618,641,4 

'  1  n  7U  1  +  1.917,800'- 


ECIJ  - 


y  /“?  (16)  (20) 


+  (21»6.g4)  (.134,090,5)  -  (.143,139,9)] 

=  .832,242,5. 


((8100)  (.100,603,3)  -  (4,900)  (.139,639,9) 


-26,244, 


(26) 


1  Cm  >■ 


■  U  ■->  •- 
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Using  (23)  E[l£]  is  next  computed 


E[IX2] 


(90)  2  +  2(200)  2  -1  200  rr 


(90)  2  +  2(100) 2  .  -1  100  f~l  (80)  2  +  2(200)  2  .  -1  200  /"7 


t  (80)  +  2(100)  *an-l  100  iZ  } 

+  - , - tan  - o-n - j 


1  r 2 


[(22,  025)  (1.262,728,8)  -  (7,025)  (1.004,044,  0) 


-  (21,600)  (1.295,153,5)  +  (6,600)  (1.055,990,3)] 


=  .149,859,2. 


Finally,  using  (24),  E[ I is  computed 


-  7  * 


l  (178)2  *  2(90) 2  ,  -1  90  rr 

l - 2 - 1-dn  "TIK 


(  178)  2  +  2(7  0) 2  „  -1  70  r^Z  (162) 2  +  2(90)2  .  -1  90  m 


,  (  162)  +  2(7  0)  -1  7  0  ^7: 

+  _ — — , - can  — j 


i  +  — ^  2  L  (11,971)  (.620,7  56,6)  -  (10,371)  (.507,553,7) 


-  (  1  0,61  1)  (.665,944,4)  +  ( 9 , 01 1)  ( .  548, 526, 7)  ] 


.693,018,1  . 


E[Pkf]  and  Var[P^f]  may  now  be  computed.  Now  using  (24)  and  (25)  one  obtains 


E[Pkf]  -  E[IX]  E[Iy] 

=  (.381,530,8)  (.832,242,5) 

=  .317,526,147  (29) 

I 

It  immediately  follows  that 

{EOkf]}  2  =  ( .  317,526,147)^ 

=  .100,822,854.  (30) 

The  expression  for  the  variance  was  obtained  as 
Var[Pkf]  *  E[ljp  E[I2]  -  (E[Pkf]}2  . 

Using  (27),  (28)  and  (30)  one  obtains 

l 

Var[Pkf]  =  (. 149,859, 2)(. 693,018,1)  -  (.100,822,854) 

=  .003,032,284. 

The  standard  deviation  of  Pkf  is 

o  =  /~003 ,032 ,284 
kf 

i 

=  .055.  (31) 

The  two-standard  deviation  confidence  interval  on  Pkf  is  immediately  obtained 
from  ( 29)  and  ( 31 )  ;  thus 

Pkf  =  .318  ±  .110  (32) 
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SECTION  IV 


ESTIMATION  OF  E[Pkf]  and  Var[Pkf]  USING 
THE  TAYLOR’S  SERIES  ESTIMATION  PROCEDURE 


1 .  Background 

Recalling  that  Pkf  is  a  function  of  the  four  input  parameters  Rx ,  R  , 

a  and  a  ,  one  may  write 
x  y 

R  R 

— -  - j -  .  (33) 

?  r~7  ? 

?o  /  R  +  ?.o 

x  y  y 

Let  R  ,  ^  ,  <7  and  o’  refer,  respectively,  to  the  mean  of  Ry ,  R  , 
x  y  x  y  y 

n  ,  and  j  .  Expanding  Puf  about  the  point  (R  ,  R  ,  o  ,  o  )  one  obtains  up  to  the 
x  y  *'  X  y  X  y 

first  order  terms: 


kf 


kf 


(MM^  ,0  )  = 


,0  )  =  P|f(R  »R  .0  ) 
kf  x’  y’  x  y  kf  x  y  x  y 


3  P 


MR  -  R 


kf 


x  "  x'  3 R 


3P 


kf 


-  -  -  -  +  (Ry"Ry^  3R 
Rx’Ry’V0y  y  y  y 


R  ,R 

x  y  x  y 


3  P 


+  j  -u 


kf 


v'x  x  ;  3o 


3P 


kf 


+  (o  -o  )  - - 

r  r  o  o  y  y  3oy 

x  ’  y  ’  x  ’  y  J 


(34) 


R  ,R  ,a  ,o 
x  ’  y  ’  x  y 


Note  that  the  partial  derivatives  of  Pkf  with  respect  to  the  four  variables  Rx , 

R  a  ,  and  a  am  to  bp  evaluated  at  the  mean  values  of  the  variables, 
y  ’  x  y 

9 .  Estimation  of  FTP  ] 

k  f 

Taking  expectations  on  both  sides  of  (34)  yields  as  a  first  approximation 


^kf'VVt 


,o  ) ]  =  P,  ,( R  ,R  ,o 
'  y’ '  kfx  x ’  y  x 
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Thus,  at  this  level  of  approximation  the  expected  value  of  is  directly 
obtained  by  setting  in  (33)  the  values  of  Rx,  Ry,  ox» and  oy  equal  to  their 
expected  value.  More  explicitly  one  can  write  using  (33) 


E[Pkf]  =  E[Pkf(Rx,Ry,ox,oy)] 


(35) 


jam?  /  IT2  +  "zo* 

x  x  y  y 


We  now  refer  to  the  definitions  of  Ix  and  Iy  as  given  in  expressions  (6)  and 
(7).  It  is  then  clear  that 


ECU  - 


I  •R2  +  Zo2 
x  x 


and 


■ 


/? 


+  2? 
y 


and  E[Pkf]  =  E[Ix]  •  E[Iy]  . 

It  should  be  noted  that  an  improved  approximation  in  the  value  of  E[Pkf]  can 
be  obtained  by  incorporating  additional  terms  in  the  Taylor's  series  expansion 
given  by  (34). 

3.  Estimation  of  Var[Pk^] 

First,  expression  (34)  is  written  as 

Pkf(Rx,Ry,Ox,Oy)  -  Pkf(^x*^y »ax»°y) 


(36) 


3P, 


kf 


-  3  Pi 


kf 


(Rx-*x)  TR — !TT  ,TT  ,o  ,o  +  (Ry"V  TO 

x'x  y  x  y  J  J  y 


1?  ,TT  ,a  ,a 
x  y  x  y 


+  (o  -a  ) 

V  y  y  / 


3P 


kf 


x~V  To“ 


3P 


kf 


Tf  ,T?  ,a  ,o  +  (°y'V  7o“ 
x  ’  y  ’  x  y  J  J  y 


1?  ,TT  ,o  ,o  ' 
x  ’  y  x  y 


Souarinn  and  taking  expectations  on  both  sides  of  (36)  yields 


3P 


Var[Pkf(Rx,Ry,o  o  )]  --  Var[R  ]  (-J^) 


x  ’  y ’  x’  y 


W 


IT  ,1T  ,o  ,o 
x  y  x  y 


j  e 

v<'V  7T 


dP 


kf 


*  Var[°x]  ^ — j 


77  ,JT  ,o  ,o 
x  y  x  y 


TT  ,JT  ,o  ,o 
x  y  x  y 


«  i  r 


:  z  :  • 

*  ✓  ’  «  / 


a?kf 

3  Cov[R  ,  R  ]  L  i-jfl—J 
7  x 


,  aPkf,  ..3Pkfn 

'°ks- 

VRy’V0y 


c  -  )  v 


d  R 


K  f 

J  **  J  J 


1  R  ,R  ,o  ,o 
x  y  x  y 


*  '  >  'XT''  V 


r  x  t 
6  1 J 


R  ,R’  ,o  ,o 
x  y  x  y 


3Pkf  9Pkf  . 

*  CovtX,aJ 

y  y  -R  , K  ,o  ,o 
x  y  x  y 


y  y  • 


3  Pl,  y  dPk  t 

2  C0Vt°x’Jy]  I'vJJ - J  (-Jo  }j 

x  y 


TT  ,1T  ,o  ,o 
x  y  x  y 


(37) 


In  the  special  case  when  Rx ,  Ry,  and  Oy  are  independently  distributed, 


expression  (37)  becomes 
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The  evaluation  of  ^ —  and  —  proceeds  similarly 

°*  °y 


3^kf  1  ?  7  ^ 

IT"  ("x>  <-2>  <4\>  (Rx+  2‘«>  -  ' 

*  /  +  2a  ^ 

y  y 

2a  R  R 

xx  y 


Hence , 


2  2  2  2 
3P.  ,  '  4a  R  R 

/  kf i  _  x  x  .  y 

3°*  "  <R2  *  *2)3  '  *  2°y>  ' 


Similarly, 


(R2*2o2)2  /RX*2»X 

y  y 


3P  f  2  4o2  R2 


1  2  2  * 

(Ry  +  &J)  RX  ' 


Expressions  (39)  through  (45)  can  then  be  substituted  either  in  (37)  or  (38) 


to  obtain  Var[Ptf(Rv,Ru,a^,aw)].  Writing  Var[Ptf]  for  Var[Pkf ( Rx ,Ry to*  ,a  ) ]  , 


v  ai  L  r  \  ^  /  J  •  " 

one  obtains,  for  example,  from  (38) 


.-4  -5-2 
4a  R 

Var[P  ]  - - — 1 - x  Var[R  1 

— o  — 2  — 2  — 2  J  * 
(R  +2a^)(RS2oM 
y  y  x  X 


. — 4  5-2 
4a  R 

+  - 1 — x -  Var[R  ] 

—7  —7  —7  —2  J  y 

(R  +2a^)(R  ,+2o„) 


(*W)(R2+Z<£)' 


4 a2  -R2  I?2 
y  x  y 


Var[oy] 


Var[Pkf]  » 


4a2  "R2 
_ x  y 


(°x2  Var[Rx]  +  Var[c¥]) 


4^2  -jj2 

- — -*-* - j  C5?  Var[R  ]  *  K*  Var[o  ]) 

( 25x)fRy+2^) 


4.  Example 


Referring  to  the  example  cited  in  Section  III-4,  the  following  values  are 
computed  for  the  necessary  data  to  compute  E[P|<f]  and  VarCP^f] 


\  -  85  ft; 


VarCR  ]  -  12*&L 


100  ft2 

“T7  rt 


■R  =  170  ft; 


Var[Ry]  ■  W*-™*:.  -  «|  ft2 


o  =  iso  ft; 


Varfo  1  =  (200-10°)2  =  1°»°00  ft2 


To  compute  ELP^f],  refer  to  expression  (35)  to  obtain 


/  (85)2  +  (2) (150) 2 


s 


(170r  +  (2)  (80)' 


E[Pkfl  -  ( 


85  ,  ,  170  s 

228. 52 J  l204.20J 


=  (.372)  (.832) 

=  .310.  (49) 

To  compute  VarfP^f]  one  can  refer  either  to  expression  (47)  or  expression  (48). 
f.xpression  (47)  is  used  since  the  contribution  of  each  variance  component  can  be 
identified. 


3 .  Contribution  of  VarCR^l 


4^  R2  Var[R  ] 


( 4 )  ( 1 50 ) 4  ( 170)2  (-i^) 


(R2  +  2a2)  (R2  +  2a2) 
x  x  y  y 

[ (85) 2  +  2(1 50 ) 2 1  [ ( 170) 2  +  2(80) 2] 

=  .000,082,1. 

Contribution  of  Var[Ry] 

4J4  R2  Var[R  ] 

\J  Y  u  V 

(4)  (80)4  ( 85) 2  (2^) 

(Rx  +  2<V(Ry  +  2aP 


[ (85) ^  +  2( 150) £][( 170;  +  2(80)d] 

.000,006,7. 


Contribution  of  Var[oxl 


4o2  R2  R2  VarTo  ] 
x  y  x  x 1 


(R2  +  2o2)  (R2  ♦  2o 2 ) 
x  x  y  y 


(4)  (  150)  2  (170)  2  ( 85) 2 


[(85)  +  2(150)1  [(170)  +  2(80)  1 

.002,636,5. 


Contribution  of  Var[o^] 


4?2  ^  nj  Var[0y]  (4)  (SO)2  (85)2  (170)2  (ig) 

- J  - - 

(R l  +  2^)(Rj  +  2^)  C(85)2  +  2(150)2][(170)2  +  2(80)2]3 

=  .000,047,1. 


The  variance  of  is  thus  given  by 

Var[Pkf]  =  .000,082,1  +  .000,006,7 

+  .002,636,5  +  .000,047,1  =  .002,772,4. 

The  standard  deviation  of  Pkf  is 

oD  =  /  .002,772,4  =  .053  (50) 

Kkf 

The  two-standard  deviation  confidence  interval  on  Pkf  is  obtained  from  (48)  and 
(49)  and  is 

Pkf  =  .310  _+  .106.  (51) 

The  results  are  summarized  as  follows 


Component  Variance  Percentage 

Range  weapon  radius  Rx  .000,082,1  2.96 

Oeflection  weapon  radius  Ry  .000,006,7  .24 

St.  Oev.  of  aiming  error  in  range  .002,636,5  95.10 

St.  Dev.  of  aiming  error  in  deflection  Oy  .000,047 ,1  1.70 


Total 


.002,772,4 


100.  00 


Comparing  the  results  of  (32)  obtained  through  the  subjective  estimation 
procedure  which  the  results  of  (50)  obtained  using  the  Taylor's  series  estimation 
procedure,  it  becomes  clear  that  the  two  methods  agree  very  closely  and  that  the 
estimates  are  robust. 

5.  Further  Considerations 

The  previous  model  could  be  adapted  to  account  for  the  following  three 
special  situations: 


x  y  x  y 


2.  Rx  t  Ry  ; 

3.  Rx  =  Ry  - 


O  =0 

x  y 

R;  o  =  a 
x  y 


a 

a  . 


The  most  interesting  of  these  cases  is  the  second  one  as  it  corresponds  to 
guided  weapons  in  which  the  aiming  error  is  the  same  in  all  directions.  It  can 
be  shown  then  that  given  that  the  weapon  is  aimed  in  the  neighborhood  of  an 
origin  point,  the  probability  that  the  weapon  will  impact  between  radius  r  and 
r  +  dr  (of  concentric  circles  centered  at  the  origin)  is  given  by  the  Rayleigh 
density  function 

1  ,  r2 

h(r)  dr  =  -y  exp[-  — ~]  r  dr  • 
a  2  a 

The  probability  of  kill  due  to  fragmentation  can  be  found  directly  from  (7) 

by  setting  a  =  a  -a.  Thus, 
x  y 


pk  f  ( Rx  ,0 ) 


/  R2  +  2a2 
x 


/7717 


(52) 


Using  the  technique  based  on  the  Taylor's  series  expansion  one  finds 


E[P„,(R„,Ru.a)]  = 


kf'  x*  y 


(5 


/  R2  +  2a2  /  R2  +  2o2 

x  y 


3P 


Var[Pkf(Rx,Ry,0)]  =  Var[Rx]  (j^-) 


R  ,R  ,o 
x*  y* 


3P 


kf, 


*  (jf1) 

J  \J  —  —  — 


R  ,R  ,a 
x*  y 


3P 


kf, 


+  Var[a]  (^— ) _ _  . 


(5 


R  ,R  ,o 
x*  y’ 


Carrying  on  the  differentiation  process,  one  obtains  similar  to  (39)  and  (41) 


3P 


kf 


2 

20  R 


3R 


1 


(5 


/  R2  +  2o 2  (R2  +  2o2)  2 
y  x 


3P 


kf 


2o2  R 


3R  _  2  ’ 

/  R2  +  2o2  (R2  +  2o 2)  2 
x  y 


To  obtain  3Pkf/do  refer  to  (52)  and  differentiate  partially  with  respect  to  o 
yield 


dp 


3 

2 


~  R  R  [.  I  (R2  +  2o2)  "  (4o)(R2  +  2o2) 

3o  xyL?'x  '  ' '  v 


3 

2 


(R l  *  2o2J  '  (-  h  (R2  +  202)  (40)] 


to 


+ 


=  -?o  Rx  Ry  (Rx  +  (Ry  +  2oc) 


!(«*  *  t°’)A  *  ("y  *  2»2)‘‘] 


-2o(R  R  )  (R2  +  R2  +  4a2) 
v  x  y'  v  x  y _ 


3  3  * 

(R?  +  2a2)  (R2  +  2a 2) 

x  y 


Substituting  (55),  (56)  and  (57)  in  (54)  yields 


4?1  R 2 


Var[Pkf]  =  - 1 -  Var[R  ] 

— ?  — ?  —2  — ?  J  * 

(R  +  ?a^)(R  +  ?„a~) 

y  x  ' 


,-A  f2 
4o  R 

+ - 5"  Var[R  1 

(R2  +  ?a2)(R2  +  2a2) 
x  y  ' 


4a2  R2  R2  (R2  +  R2  +  4o2) 

♦ - - * - L  Var[a] 


(R2  +  2a2)  ( R2  +  2a2) 


Lxampl i 


In  example  III -4  assume  the  weapon  to  be  guided  and  let  the  parameters  of 


interest  be:  Rx  =  85  5  ft,  R  =  170  +_  8  ft  and  a  =  30  _+  5  ft.  It  is  required 


to  determine  Pl 


Rx  =  85  ft,  Var[Rxl 


100  f  2 
12 


R  =  170  ft  ,  VarrR  ]  = 

y  y 


256  2 

12  1 


a  =  30  ft ,  Var[o] 


100  f  2 
12  TZ 


Using  (53)  yields 


E[Pkf]  -  - . -  f—  .  •  - .  *70-- 

J  (85)2  +  2(30)2  /  (170)2  +  2(30)2 

•  (.895) (.970)  =  .868  .  (59) 

Expression  (58)  is  now  used  to  determine  the  contribution  of  the  various  variance 
components  to  the  total  variance  which  is  Var[Pkf]. 

a.  Contribution  of  Var[Rxl 

4^  R2  Var[Rxl  (4)  (30)4  (170)2 

__  _  3  = - 3 - 

(R2  +  2a2)  (R2  +  2a2)  [ ( 35 ) 2  +  (2)(30)2]  [(170)2  +  (2)(30)2] 

=  .000,034,6. 

b.  Contribution  of  Var[R  ] 

4^  R2  Var[Ry]  (4) (30)4  (85)2  (^||) 

(R?  +  2o2)(R2  +  2o2)3  [ ( 85 ) 2  +  (2)(30)2][(170)2  +  (2)(30)2) 

*  y 

=  .000,001,9. 


c.  Contribution  of  Var[a] 


4a2  R2  R2  (R2  »  R2  ♦  4a2)  Var[g] 
(*l  +  2a2)  (R2  +  2a2) 


(4)(30)2  (85)2  ( 1 70 ) 2  r ( 85 ) 2  +  (170)2  +  (4)(30)2]'  (^|) 

3  1  3 

t(85)2  +  (2)(30)2]  [( 170)2  +  (2) ( 30)2] 


=  .000,464,8. 
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The  variance  of  Pkf  is  thus  given  by 


Var[Pkf]  =  .000,034,6  +  .000,001,9  +  .000,464,8 
=  .000,501,3. 

The  standard  deviation  of  Pkf  is 


o  =  /  .000,501,3  =  .022.  (60) 

Pkf 


The  two-standard  deviation  confidence  interval  on  Pkf  is  obtained  from  (59)  and 
(60)  and  is 

Pkf  =  .868  +_  .044. 

The  variance  results  are  summarized  as  follows: 


Component 

Range  weapon  radius  Rx 
Deflection  weapon  radius  Ry 
St.  Dev.  of  aiming  error  a 


Variance 

000,034,6 

000,001,9 


Percentage 

6.90 

.38 

92.72 


Total 


.000,464,8 

.000,501,3 


100.00 


SECTION  V 


CONCLUSIONS 

In  this  report,  two  methods  are  developed  to  measure  the  variability  of  the 
probability  of  kill  P^f  due  to  fragmentation  of  weapons  delivered  from  air.  It 
is  assumed  that  delivery  error  is  present  in  killing  the  target,  both  weapon  and 
target  being  idealized  as  points.  A  two-parameter  Carleton  damage  function  is 
used  to  describe  mathematically  as  a  function  of  the  distance  between  the 
location  of  the  target  and  the  location  of  weapon  explosion. 

It  is  shown  that  it  is  possible  to  provide  two-standard  deviation  confidence 
intervals  on  Pkf  and  that  the  estimates  are  robust, yielding  approximately  the 
same  result  for  the  two  methods. 

It  is  suggested  that  future  work  investigates  the  case  of  the  more  general 
three-parameter  Carleton  damage  function  using  the  same  methodology  to  compute 
confidence  intervals  on  P^. 
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